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Abstract 

We study to which extent SUSY extensions of the Standard Model can describe 
the excess of events of 3.0 standard deviations observed by ATLAS in the on-Z signal 
region, respecting constraints by CMS on similar signal channels as well as constraints 
from searches for jets and GMSB-like scenarios are typically in conflict with 

these constraints, and do not reproduce well the shape of the distribution of 

the data. An alternative scenario with two massive neutralinos can improve fits to the 
total number of events as well as to the and distributions. Such a scenario 

can be realised within the NMSSM. 


1 Introduction 


After the first run of the LHC at a center of mass (c.m.) energy of mostly 8 TeV, no 
signihcant excesses have been observed in searches for physics beyond the Standard Model 
HE]. These searches cover a wide range of possible signatures, notably various combinations 
of jets, missing transverse energy and leptons (electrons or muons). 

Same-flavour opposite-sign dileptons can be classihed into “off-Z” leptons (typically with 
an invariant mass mu < 81 GeV or mu > 101 GeV), and “on-Z” leptons with 81 GeV < 
mu < 101 GeV. Often, leptons and in particular on-Z dileptons are vetoed in order to 
suppress Standard Model (SM) backgrounds. On the other hand, some decay cascades 
of supersymmetric (SUSY) particles could be particularly rich in off-Z dileptons (in the 
presence of light sleptons), or on-Z dileptons if Z bosons appear particularly frequently in 
these cascades. 

Recently, results of searches for SUSY particles in events with dileptons, jets and 
have been published by the GMS and ATLAS collaborations The aim was to test 

scenarios of gluino pair production in which the gluinos g decay via sleptons (leading to 
off-Z dileptons), and scenarios of gauge mediated SUSY breaking (GMSB) or generalised 
gauge mediation (GGM) where the gluinos undergo 3-body decays into quark pairs and 
a neutralino Xv The latter may decay subsequently into a nearly massless gravitino G 
and a Z boson, leading to on-Z dileptons. The corresponding gluino decay chain is then 
g —)■ q + q + Xi q + q + Z + G. Relevant parameters are the gluino mass m^, the neutralino 

mass m^o, and the branching fractions of the involved decays. 

Whereas no significant excesses were observed by GMS in [3] (up to an excess of 2.6 stan¬ 
dard deviations in the dilepton mass window 20 GeV < mu < 70 GeV), an excess of 
3.0 standard deviations was reported by ATLAS in |1] in the on-Z signal region; Summing 
electron and muon pairs, 29 events passing the cuts were observed versus 10.6 ± 3.2 back¬ 
ground events expected. No attempt was made in |1] to explain the excess in terms of a 
specific model; instead, weaker exclusion limits than expected were shown in the m^ — m^o 
plane of GGM models. Various studies of scenarios which could contribute to this excess 
have recently been published [511T3]. 

Z bosons decay dominantly hadronically. Thus, whenever gluinos are pair produced, 
in most cases each of the two gluino cascades will produce no dileptons, but two hard 
jets: either from q + q ii m^ ^ m^o > Mz, or from hadronic Z decays if m^ > m^o ^ Mz 
implying a neutralino much heavier than the gravitino, i.e. energetic Z bosons. Hence, both 
scenarios are subject to constraints from “standard” searches for SUSY in events with hard 
jets and [miTS] . even if one considers simplified models where squarks are assumed 

to be decoupled and gluino pair production is the only process taken into account. 

In order to study the impact of these constraints on GMSB-like scenarios, we simulated 
various conhgurations of gluino and Xi masses. Using the latest version 1.2.0 of GheckMATE 
[16] we found that constraints from HU (a preliminary version of mi) on final states with 
jets and are very restrictive, and supersede even the recent GMS constraints from [3] 

in the m^ — m^o plane. Exceptions are scenarios with reduced branching fractions for the 
considered decay chain, without allowing for alternative final states leading to jets and 

Tpmiss 

H/'-£ . 

In the present paper we study to which extent a scenario with two heavy neutralinos in 
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the gluino decay cascade can contribute to the ATLAS signal region, circumventing con¬ 
straints from searches for jets and The gluino decay cascade considered subsequently 

is of the form 


g^q + q + xl^Q + Q + Z + xl 


( 1 . 1 ) 


with 


and sketched in Fig. 1. 


^ m-g, 


m^o ~ m^o — 100 GeV 

A1 A2 


( 1 . 2 ) 



Figure 1: Gluino decay cascades involving two neutralinos X 2 Xi- 

Now jets from both steps of the gluino decay cascade (including the jets from the 
Z boson) are relatively soft, and constraints from searches for jets and are easier to 

satisfy unless the mass splitting rrig — m^o is too large. Such a scenario has been considered 
recently also in [9]. We will compare their results to ours in the conclusions. 

In the following we consider hrst simplihed models with 100% branching fractions for 
both steps of the gluino decay cascade. We simulated corresponding events, verihed which 
scenarios satisfy the constraints from the GMS [3] and other SUSY searches, and applied 
the cuts of ATLAS |1]. We will compare the signal rates and various distributions to the 
data given in |1] , and to a GMSB-like scenario (the latter with reduced branching fractions 
in order to comply with constraints). Gonstraints from GMS [3] prevent an excess as large 
as 3.0 standard deviations in the ATLAS signal region, but about 14 signal events on top 
of the background are possible. 

However, the question arises in which SUSY scenario such a neutralino spectrum and, 
notably, such a dominant gluino decay cascade are possible: What can prevent a dominant 
g —)■ q + q + Xi decay which is favored by phase space? In GMSB the role of Xi is played by 
the nearly massless gravitino, which has tiny couplings to the MSSM-like sparticles and is 
not produced unless, due to R-parity conservation, it is the only decay channel. A heavier 
neutralino Xi with small couplings to the MSSM-like sparticles, as required in the present 
scenario, is possible in the NMSSM [18] in the form of the singlino, the fermionic partner 
of the singlet superheld S whose vacuum expectation value generates dynamically a p-term 
(a SUSY mass term for the two Higgs doublets in the MSSM) of the order of the SUSY 


2 




breaking scale. We find that there exist indeed scenarios within the parameter space of the 
NMSSM for which the gluino decay cascade in Eq. fll.ip is dominant. 

In the next section we describe details of the simulation and cuts. Results for simplified 
models and the description of a NMSSM scenario are given in section 3. We conclude in 
section 4. 


2 Simulations and cuts 

We have simulated events at the LHC at 8 TeV using MadGraph/MadEvent [I9] which in¬ 
cludes Pythia 6.4 [20] for showering and hadronisation. The emission of one additional hard 
jet was allowed in the simulation in order to obtain realistic distributions for kinematical 
variables. The production cross sections were obtained by Prospino at NLO [2IIE2]. 

First, the output was given to CheckMATE version 1.2.0 [T6| which includes the detector 
simulation DELPHES [23] and compares the signal rates to constraints from various search 
channels of ATLAS and CMS. All searches present in CheckMATE version 1.2.0 have been 
verified; the most relevant ones (with the largest ratio for the event yield to where 
is the observed 95% CL upper bound) are obtained from the ATLAS search for jets and 
in [H]. 

Second, the Pythia output was given directly to DELPHES and analysed according to 
the object identification and selection criteria given in [3lll], respectively, and finally the 
corresponding cuts were applied. 

For the ATLAS on-Z searches [1] these were as follows: > 225 GeV; > 2 jets with 

Pt > 35 GeV; two same-flavour opposite-sign leptons with px > 25 GeV for the leading, 
Pt > 10 GeV for the sub-leading lepton; > 600 GeV where iLx = -|- 

(including jets with px > 35 GeV); and finally 81 GeV < mu < 101 GeV. 29 events 
passing the cuts were observed, whereas 10.6 ± 3.2 background events were expected. 

For the events passing the cuts, distributions of mu, and the jet multiplicity 

Ajets were shown in [1] separately for the electron and muon channels. These distributions 
were compared with those expected from two GGM benchmark points with gluino masses 
and neutralino masses of {mg,m^o) = (700, 200) GeV, (900, 600) GeV, respectively. We 
found, however, that both points violate constraints from ini on final states with jets and 

jpmiss 

We compared the expected properties of the two GGM benchmark points in [1] to the 
results of our simulation and found that they agree within ~ 30%. We conclude that 
the results of our simulations deviate by a systematic error of up to ~ 30% from the 
more realistic (detector-) simulation of the experimental collaboration. We can expect that 
this systematic error cancels to a large extent when comparing the properties of different 
simulated scenarios, but should be taken into account when comparing to the actual data 
from [1]. Since it is of the same order (actually somewhat larger) than the difference in 
the acceptances of dielectrons and dimuons in [1], we found it reasonable to consider the 
sum of the data of dielectron and dimuon events not only for the signal rate, but also for 
the kinematical distributions and the expected SM background in order to obtain a larger 
statistics. 

In the GMS on-Z searches [3], no cuts on iLx were applied. Signal jets were required 
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to have px > 40 GeV. Six and A^jets-dependent on-Z signal regions were defined: 

£rniss _ ^QQ _ 200, 200 — 300, > 300 GeV and iVjets >2, >3, respectively. Finally the 
GMS and ATLAS analyses differ slightly in the jet algorithms and in the lepton acceptances. 
Gomparing the signal rates obtained by our simulations of the two GMSB-like benchmark 
points to the simulations in [3] we found again that they agree within ~ 30%. 

As already stated above, no significant excesses were observed in the on-Z signal regions 
by GMS. Hence the event yields in the six on-Z signal regions lead to constraints on any 
scenarios which attempt to explain the ATLAS excess. In the next section we discuss 
by means of benchmark points to which extent the ATLAS excess can be matched in 
consideration of these constraints, as well as constraints from HU on final states with jets 
and 


3 Results 

First we considered GMSB-like simplified models with a branching fraction of 100% for the 
g^q + q + ')(^^q + q + Z + G decay chain. Then, however, constraints from the search 
for jets and in [T7] as tested by GheckMATE (TB] require rrig > 1050 GeV for small 

m^o ~ 150 GeV, and larger gluino masses for larger m^o. Accordingly contributions to the 
ATLAS signal region cannot exceed ~ 5 events. Moreover the distribution of and 

Ht peak towards large values (most events have Ht > 1500 GeV) in sharp contrast to the 
data in |1]. 

In realistic models, the branching fractions for the steps of the above gluino decay chain 
can well be below 100%. Below we will consider a GMSB-like benchmark point “GMSB” 
with {rrig, m^o) = (800, 600) GeV and a branching fraction of 10% for the above decay chain; 
such a small branching fraction makes it compatible with the GMS constraints. (A heavy 
Xi was chosen in order to shift the peak of the distribution towards lower values.) For 
the remaining 90% of the gluino decays one has to expect that, depending on the complete 
spectrum and branching fractions, they contribute to the signal regions in the search for 
jets and in [1^. One can make the somewhat optimistic assumption that these 

contributions do not exceed 50% of the contributions of the g ^ q + q + Xi q + q + Z + G 
decay chain. Then this point remains within the constraints from HZ], but contributes 
about 10 events to the ATLAS on-Z signal region. 

Next we consider simplified models with two heavy neutralinos whose decay chain is 
depicted in Fig. 1. Assuming a branching fraction of 100% for this decay chain, gluinos can 
be as light as 800 GeV without conflict with constraints from the search for jets and in 
m - under the condition, however, that m^o and m^o are relatively large such that all jets 
remain relatively soft. We studied two benchmark points PI and P2 with {rrig, m^o,my^o) = 
(800,790,690) GeV and (m^, m^o, m^o) = (800,600,500) GeV, respectively. Such scenarios 
belong to the few exceptions allowing for gluinos with a mass below 1 TeV, see the study 
in [21]. 

For PI with nig — m^o = 10 GeV the jets from the first step g ^ q + q + X2 the 
decay cascade are very soft, as are the jets from Z decays from the second step X2 
Z + Xv Practically all energy of a single gluino decay cascade goes into However, 

for typical kinematical configurations the momenta of Xi tend to be back-to-back in the 
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transverse plane, leading to a reduction of of the complete event. Only for relatively 

rare kinematical configurations (and/or extra jets from initial state radiation as included 
in our simulation), of the complete event can assume large values. For P2 with 

rrig — m^o = 200 GeV the jets from the first step g ^ q + q + X2 ^he decay cascade 
are harder, leading to less One aim is to study the impact of this difference on the 

distributions of kinematical variables. 

For all benchmark points we assumed practically decoupled squarks with masses of 
3 TeV; then the gluino pair production cross section from prospino at NLO is 128 fb. (Since 
stops and sbottoms are assumed to have masses of 3 TeV as well their pair production 
does not contribute to the signal.) We deliberately chose identical gluino masses for all 
points in order to maintain a common production cross section; therefore all differences in 
contributions to signal regions and kinematical distributions originate from the neutralino 
sector. The masses of the latter are recalled in Table 1 below. 

In addition we indicate in the Table 1 in how far the benchmark points GMSB, PI 
and P2 satisfy constraints from the six signal regions of the GMS on-Z searches in [3] 
(including 30% systematic errors from the simulation). The 95% GL upper limits for the 
six signal regions of the GMS on-Z searches had already been obtained in [9]. We hnd that 
the central values of event yields of the benchmark points are below these 95% GL upper 




GMSB 

PI 

P2 

Gluino/neutralino masses 





m~g 


800 

800 

800 

m^o (GMSB), m^o (PI, P2) 


600 

790 

600 

rua (GMSB), m^o (PI, P2) 


0 

690 

500 

Gonstraining signal regions 

C95 

^obs 




GMS, N^ets > 2, 100 < < 200 

207 

2.0 ± 0.6 

19.4 ± 5.8 

55.7 ± 16.7 

GMS, N^ets > 2, 200 < < 300 

20 

2.6 ± 0.78 

8.1 ± 2.4 

23.7 ± 7.1 

GMS, Njets > 2, 300 < 

7.6 

7.0 ± 2.1 

6.1 ± 1.8 

7.0 ± 2.1 

GMS, N^ets > 3, 100 < < 200 

89 

1.9 ± 0.57 

8.5 ± 2.6 

48.4 ± 14.5 

GMS, N^ets > 3, 200 < E™*" < 300 

16.1 

2.4 ± 0.72 

4.7 ± 1.4 

21.1 ± 6.3 

GMS, Njets > 3, 300 < E™"" 

8 

6.4 ± 1.9 

3.8 ± 1.1 

6.7 ± 2.0 

ATLAS, GT 

2.4 

0.73 ± 0.22 

0.79 ± 0.24 

1.31 ± 0.39 

ATLAS, EM 

28.6 

21.7 ± 6.5 

1.32 ± 0.39 

15.8 ± 4.7 

ATLAS, ET 

8.3 

7.8 ± 2.3 

0.70 ± 0.21 

4.13 ± 1.24 

ATLAS on-Z SR (obs. excess 18.4) 


9.8 ± 2.9 

5.6 ± 1.7 

13.6 ± 4.1 


Table 1: Sparticle masses of the benchmark points GMSB, PI and P2 (in GeV), event yields 
including 30% systematic errors from the simulation of the benchmark points GMSB, PI 
and P2 in the six signal regions of the GMS on-Z searches in [5], and in the most constraining 
signal regions GT, EM and ET of the ATLAS search [1^. The ranges of for the six 
GMS signal regions are given in GeV. The last line indicates the contributions to the ATLAS 
on-Z signal region. 
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limits with the exception of P2 in the bins N-^^ts ^ 2, 200 < < 300 and A^jets > 3, 

200 < < 300. However, taking the systematic errors from the simulation into account, 

the CLg = CLg+b/CLh values for P2 in these bins are 0.11 and 0.09, respectively, i.e. well 
above the 95% CL exclusion limit of 0.05. 

Out of the 10 signal regions in the ATLAS search m for jets and we show 

the event yields for the signal regions CT, EM and ET which give the largest ratio event 
yield/S'^^g for the points PI, P2 and GMSB, respectively. All these signal regions require 
^ QgY^ ^ Q^y fgj, leading jet, and pt > 60 GeV for 3 additional jets 
(CT), Pt > 60 GeV for 5 additional jets (EM and ET). EM and ET differ by mes{incl.) > 
1200/1500 GeV, respectively (see [T7] for more details). 

We recall that the event yields for the point GMSB assume only a branching fraction 
of 10% into the considered gluino decay chain. In Table 1, only the contributions from the 
simulated decay chain are shown. Within the systematic error bars, 50% more events from 
other gluino decays are allowed to contribute to the signal regions in the ATLAS search Id 
for jets and in order to saturate the bound from the signal region ET. 

Finally we compare the contributions of the benchmark points GMSB, PI and P2 to the 
ATLAS on-Z signal region, summing dielectrons and dimuons, in the last line of Table 1. 

We see that a price has to be paid for the very compressed gluino—y® ~ Xi spectrum in 
PI: Due to the softness of the jets, not enough jets satisfy the cut Ajets > 2. The GMSB 
point seems to do quite well, despite its gluino branching fraction being reduced by a factor 
~ 1/10. The best ht is given by P2 with its less compressed gluino—X 2 ~ Xi spectrum. 

Next we consider the distributions of kinematical variables. As stated above we combine 
the ATLAS dielectron and dimuon data (despite the different acceptances) in order to 
enhance the visibility of possible trends. We only show the (dominant) statistical error of 
the data; we are not in a position to combine the partially correlated systematic errors. 
In the figures below we show the data with the expected SM background contribution 
subtracted, with the aim to expose possible desirable features of signal contributions (see |1] 
for the error attributed to the expected background). 

We start with in Fig. 2 where we compare the data with the expected background 

subtracted to the GMSB scenario and with the two heavy-neutralino benchmark points 
PI and P2. We simulated 500.000 events for each scenario. Each expected event for the 
LHG run I as shown in Fig. [2] corresponds to 10 simulated events, which allows to estimate 
the statistical errors. These are smaller than the estimated systematic errors from our 
simulation, and much smaller than the statistical error of the data. 

The measured event numbers seem to decrease continuously with E™*®® (within the 
error bars, and note that the rightmost bin includes the overflow), whereas the E^®® 
distributions of the GMSB and PI points are nearly flat: In these scenarios, nearly all 
energy is transformed into missing energy which prefers accordingly large values of Ej?*®®. 
(Note that E™®® is shown after the application of all cuts, notably on Ht > 600 GeV. For 
PI with its compressed spectrum this cut selects atypical kinematical configurations with 
particularly large E™®®.) 

For a quantitative comparison we compute the reduced statistic 

, 1 - Ns(i)f 

iVbin. - 1 ^ ' 
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Figure 2: Comparison of from the data in [1] (with the expected background sub¬ 

tracted) to the benchmark points GMSB, PI and P2 dehned in the text. Error bars on the 
data are statistical only. The rightmost bin includes the overflow. 


for each benchmark point, where Nd-b{i) is the data with the expected background sub¬ 
tracted (as shown in Fig. 2). a{i) combines the statistical error of the data shown in Fig. 2 
and the systematic error of 30% of our simulation (with respect to which the systematic 
error of the background is negligible). 

We obtained Xred = 0-69 for GMSB, Xred = PI Xred = 0-61 for P2. Hence 

the scenario P2 with its larger splitting between the gluino and the xt rnasses describes 
best the shape of the distribution. Of course, the scenario P2 prohts also from its 

larger total event rate. 

In Fig. [3] we compare the data on Ht (with the expected background subtracted) with 
the GMSB scenario and the benchmark points PI and P2. Since represents most of the 
visible transverse energy, the point PI with its compressed spectrum peaks at low values 
of Ht- This coincides with the trend of the data, but the total signal rate (limited by 
constraints from GMS) is small, as indicated in Table 1. 

For the reduced statistic we find Xred — for GMSB, Xred — 9-69 for PI and 
Xred = 0-36 for P2. Again, the benchmark point P2 provides the best agreement with the 
shape of the distribution despite its somewhat less compressed spectrum. 

Finally we turn to the distribution of the jet multiplicity in Fig. 01 The trend of the data 
towards low jet multiplicities is reproduced only by PI with its excessively low signal rate. 
The jet multiplicity of simulations is sensitive, amongst others, to the matching between 
soft and hard QGD radiation, accordingly this quantity has to be considered with some 
reserve. 

For the reduced statistic we find y^g^ = 1.03 for GMSB, y^g^ = 1.08 for PI and 
Xred = 1-57 for P2. In this case the trend of the data is not well reproduced by the point 
P2. But since the scenario P2 provides the best fit to the ATLAS signal rate and the 
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Figure 3: Comparison of from the data in |1] (with the expected background subtracted) 
to the benchmark points GMSB, PI and P2 dehned in the text. Error bars on the data are 
statistical only. The rightmost bin includes the overflow. 
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Figure 4: Comparison of the jet multiplicity from the data in |1] (with the expected back¬ 
ground subtracted) to the benchmark points GMSB, PI and P2 dehned in the text. Error 
bars on the data are statistical only. 

distributions, it would be interesting to know about SUSY extensions of 
the Standard Model which share the features of this simplihed model. As discussed in the 
introduction, this is possible within the NMSSM. 

Using the spectrum generator NMSSMTools [2511^ with decay branching fractions 
computed by NMSDECAY [2^ (based on HDECAY j2H]) we found that the following 















region of the parameter space of the Za-invariant NMSSM shares the following properties 
with the point P2: 

• Heavy (decoupled) squarks in order to satisfy constraints from searches for events 
with jets and in the presence of a gluino with a mass of 800 GeV. 

• A bino-like neutralino X 2 with a mass of 600 GeV, but winos and higgsinos slightly 
heavier than the gluino. (The running gaugino masses do not satisfy SU(5)-like rela¬ 
tions at the GUT scale.) 

• A singlet-like neutralino Xi with a mass of 500 GeV. Then the branching fraction for 
the decay X 2 -|- Z is 100%. 

• The loop-induced gluino two-body decay g g + Xi should be suppressed, since it 
would not contribute to the signal. It is induced by the higgsino component of 
and can be of similar order of the desired gluino three-body decay g —)■ q + q + X 2 - 
The singlino-higgsino mixing is proportional to the NMSSM-specific Yukawa coupling 
A [18], and A should not exceed ~ 0.3. (The loop-induced gluino two-body decay 
g ^ g + X2 leads to similar signals a.s g ^ q + q + X2 j but it can be expected that it 
would improve the jet multiplicity distribution shifting it towards smaller values.) 

The remaining parameters can be chosen to obtain a Standard Model-like Higgs boson 
with a mass of ~ 125 GeV. We have checked that a corresponding point in the parameter 
space with all squark masses of 2.5 TeV (leading to a gluino pair production cross section 
of ~ 150 fb in order to compensate for a gluino BR into q + q + X 2 slightly below 100%), 

tan/3 = 3.75, peff = 800 GeV, A ~ 0.28, k ~ 0.087, Aa ~ 2.7 TeV and -50 GeV 

(see [18] for the dehnitions of the latter parameters) has the properties of P2 and would not 
be distinguishable from P2 regarding the different observables shown above in Figs. [2H11 


4 Summary and conclusions 

We studied to which extent SUSY extensions of the SM can describe the excess of events 
observed by ATLAS in the on-Z signal region, respecting constraints by GMS on similar 
signal channels as well as constraints from searches for jets and For viable scenarios 

we compared the distribution of kinematical variables to the data, combining dielectron 
and dimuon events. 

Due to hadronic Z-decays, GMSB-like scenarios are typically in conflict with constraints 
from searches for jets and Assuming a 100% branching fraction for the gluino decay 

cascade including the Xi G-l-Z decay, these scenarios become viable only if the gluino has 
a mass above ~ 1.05 TeV implying a small contribution (< 5 events) to the ATLAS signal 
region. Reducing the branching fraction (including the x? G-l-Z decay) to ~ 10%, lighter 
gluinos with rrig ~ 800 GeV may contribute signihcantly to the signal region, remaining 
within the 95% GL limit of GMS. However, the Ht and notably the distributions do 
not coincide well with the trends of the data. Therefore we studied alternative scenarios 
with two massive neutralinos X 21 Xi- In order to compare the impact of different neutralino 
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spectra to GMSB-like scenarios and among themselves for fixed gluino pair production cross 
sections and gluino masses, we fixed the latter also to 800 GeV. 

A very compressed g — X 2 ~Xi spectrum reproduces somewhat better the trend of the 
distribution, but does not improve the shape of the distribution. In particular, the 

contribution to the ATLAS signal region cannot be enhanced signihcantly while remaining 
within the 95% GL limit of GMS. 

A less compressed g — X 2 ~Xi spectrum provides the best £t to the total number of events 
in the ATLAS in the on-Z signal region, as well as to the Lfx and distributions; only 

the jet multiplicity is still not well reproduced. (Larger g — X 2 mass splittings as assumed 
here become again sensitive to constraints from searches for jets and We found that 

such a scenario can be realised within the NMSSM. 

A somewhat different approach has recently been persued in [9], where the space of 
the two lightest neutralino masses within the NMSSM was scanned systematically in order 
to maximise the contribution to the ATLAS on-Z signal region respecting existing con¬ 
straints. The shapes of the kinematical variables have not been studied, however. Still, 
their main results coincide with ours: Whereas compressed spectra make it easier to satisfy 
constraints from other SUSY searches, the contributions to the ATLAS on-Z signal region 
are suppressed as well. For gluino masses below ~ 800 GeV, only a small corner in the 
plane of the two lightest neutralino masses survives the 95% GL limits of GMS. Within this 
corner (for g, Xi masses of 650, 565 and 465 GeV, respectively) the authors found a 
maximal contribution of about 11 events to the ATLAS on-Z signal region. 

Moreover, the authors of [9] considered constraints from signal regions in [T3] (2jW and 
4jW) which are not implemented in the GheckMATE version 1.2.0 [IB] used here. The 
authors applied these constraints to our benchmark point P2 and obtained a ratio for the 
event yield/F^^g of 1.19, i.e. about 20% too large, but within the systematic errors from the 
simulation. A similar excess holds for this point actually also for two GMS signal regions 
considered in Table 1. We recall that identical gluino masses of 800 GeV were chosen for all 
points to simplify comparisons. A slightly heavier gluino mass of ~ 825 GeV would reduce 
the gluino pair production cross section by ~ 20%, but with little changes in the decay 
kinematics if all mass splittings remain the same. Then this modihed point would pass all 
constraints without the help of systematic error bars, but its contribution to the ATLAS 
on-Z signal region would drop to ~ 11 events. This number coincides with the maximum 
found in [9] for the slightly different point above. 

Glearly, if the excess observed by ATLAS indicates the presence of particles beyond the 
Standard Model, it should become more visible in both ATLAS and GMS experiments at 
the run 11 of the LHG. But since it is present in the ATLAS analysis of the available data 
from run 1 we found it appropriate to discuss possible interpretations. 

Within the class of models considered here, hts to the event numbers and shapes of the 
ATLAS on-Z can be improved with respect to the GMSB scenarios considered in [1] . How¬ 
ever, perfect hts would lead to unacceptable tensions with constraints from other searches. 
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